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r-| ■ Abstract 

The isospin violating dark matter (IVDM) scenario offers an interesting possibility to reconcile 
conflicting results among direct dark matter search experiments for a mass range around 10 GeV. 
r^ ■ We consider two simple renormalizable IVDM models with a complex scalar dark matter and a 

o 

CO 



Dirac fermion dark matter, respectively, whose stability is ensured by the conservation of "dark 
matter number" . Although both models successfully work as the IVDM scenario with destructive 
interference between effective couplings to proton and neutron, the dark matter annihilation cross 
section is found to exceed the cosmological/astrophysical upper bounds. Then, we propose a simple 



scenario to reconcile the IVDM scenario with the cosmological/astrophysical bounds, namely, the 
IVDM being asymmetric. Assuming a suitable amount of dark matter asymmetry has been gen- 
erated in the early Universe, the annihilation cross section beyond the cosmological/astrophysical 
upper bound nicely works to dramatically reduce the anti-dark matter relic density and as a result, 
the constraints from dark matter indirect searches are avoided. We also discuss collider experi- 
mental constraints on the models and an implication to Higgs boson physics. 
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I. INTRODUCTION 



Light weakly interacting massive particles (WIMPs) with a mass around 10 GeV has been 
currently a subject of interest, motivated by some recent results in direct dark matter (DM) 
detection experiments. DAMA/LIBRA has claimed detections of the annual modulation 
signal by WIMPs |l|. CoGeNT has found an irreducible excess [2| and annual modulation |3j. 
CRESST has observed many events that expected backgrounds are not enough to account 
for |4| . However, these observations are challenged to the null results obtained by other 
experimental collaborations, CDMS [5(, XENON10 |(|, XENON100 p, b| and SIMPLE |g|. 

Light WIMPs have been investigated for a dark matter interpretation of those data. For 
instance, very light neutralino in the Minimal Supersymmetric Standard Model (MSSM' 



10. 
11| and the Next-to-MSSM (NMSSM) [l2|, [l3| or very light right-handed sneutrino [bj, [lj] 
in the NMSSM. On the other hand, the Fermi-LAT collaboration has derived constraints on 



s-wave annihilation cross section of a WIMP based on the analysis of gamma ray flux 16]. 
Annihilation modes of a light WIMP is now severely constrained. 

The isospin violating dark matter (IVDM) [17| has been proposed as a way to reconcile 
the tension between inconsistent results among the direct DM detection experiments, since 
different nuclei for target material have been used in the detector of each experiments. The 
possible consistency between DAMA, CoGeNT [2J and XENON jd, Lzf| was pointed out |l7| . 
while the discrepancy between CoGeNT and CDMS can not be resolved by IVDM because 
both of them use germanium as the target. However, recently it is reported 18j that CDMS- 



II Si have observed three events and its possible signal region overlaps with the possible 



CoGeNT signal region analyzed 
examined by several groups 19}- 
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ig|. The fitting data with IVDM have been 



25l |. and constraints from indirect 26] and direct [27j DM 



detection experiments also have been derived. 

In this paper, we consider two simple IVDM models with a complex scalar DM and a 
Dirac fermion DM, respectively. In most of previous works, the IVDM models have been 
proposed by introducing a new U(l) gauge symmetry with Z' boson 28|, |29( or an extension 
of the Higgs sector 29M31| . In contrast to those models, to realize the different cross section 
with respect to up quark and down quark, we introduce fourth generation quarks in the 
scalar DM model and scalar quarks in the fermion DM model, respectively. Our models are 
similar to a model briefly mentioned in Ref. [17]. 



TABLE I: Particle contents for the model S 



Fields 


SU(3) C 


SU(2) L 


U(l)y 


U(1)g 
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+2/3 
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-1/3 


+1 


4> 
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1 
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1 
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The paper is organized as follows. In the next section, we describe our models of scalar 
and fermion DMs. In section III, we identify the allowed region of the mass and couplings 
of the mediator quarks or scalar quarks by imposing the condition of the isospin violating 
elastic scattering cross section with nuclei. In section IV, we calculate the annihilation cross 
section of the IVDM to examine the resultant thermal relic density as well as the constraint 
from Fermi-LAT data for the parameter region found in section III. Constraints from collider 
experiments are discussed in section V. Section VI is devoted to conclusions. 

II. MODELS 

A. Model of scalar dark matter with fermion mediators (Model S) 



First, we consider a simple model with a complex scalar dark matter, whose particle 
contents are given in Table HI In addition to the Standard Model (SM) particle contents, we 
have introduced the SM SU(2) singlet Dirac fermions (U and D) whose representations are 
the same as SU(2) singlet up and down quarks, a complex scalar DM (</>), and a real scalar S, 
with a global U(1)g symmetry. The stability of <fi is ensured by the global U(1)g symmetry 
assumed to be conserved. All the SM particles are neutral under the global symmetry. 

The gauge and global symmetric Lagrangian relevant to our discussion is given by 



£ D -M V UU - M D DD - (fuU L </m R + f D D L $d R + h.c.) - V(H, 0, S), 



;i) 



where H is the SM Higgs doublet, u R (cLr) is the SM right-handed up (down) quark singlet, 
and V is a scalar potential for H, and S. 

We assume a suitable scalar potential for our discussion: not only the Higgs doublet but 



also the scalar S develop vacuum expectation values and we expand these scalar fields as 

H = ( ° | , S = v s + s, (2) 

with the vacuum expectation values, v = 246 GeV and v s . After the electroweak symmetry 
breaking, the SM singlet scalar and the Higgs boson have a mass mixing such that 

cos a sin a \ ( hi \ 

U ' (3) 

- sin a cos a I \ hi J 

where h\ and h% are the mass eigenstates with masses m^ < m^, respectively. The exis- 
tence of a li ght scalar particle mixed with the SM Higgs boson is constrained by the LEP 




experiments |32l. l33l] . We consider a small mixing, for example, sin a < 0.1, so that the mass 
eigenstate h% (h?) is almost the SM singlet scalar (the SM Higgs boson). For such a small 
mixing, the lower mass bound on hi disappears, and in the following analysis we consider 
m/ il < 10 GeV. We also assume that the scalar potential includes a triple scalar coupling 
parameterized as 

V D Xvscfrcf) = \v(cosahi + smah2)(j> (j>, (4) 

with a dimension- less coupling A. 

B. Model of fermion dark matter with scalar mediators (Model F) 

Next, we consider a simple model with a Dirac fermion DM, whose particle contents are 
given in Table [TTl In addition to the SM particle contents, we introduce color triplet scalars 
(Ql, Ur and Dr) which are analogous to the scalar quarks in the MSSM, and a Dirac 
fermion DM (ip). Similarly to the Model S, a global U(1)q symmetry has been introduced 
to ensure the stability of the Dirac fermion DM. All the SM fields are neutral under the 
global symmetry. 

The relevant part of Lagrangian is given by 

C D -m^ - M%Q\Q L - M*U R U R - M 2 D D j R D R 
+AuQ{HU R + A d qIhD r + h.c. 
-MQlQl - faJptJ^UR - f Rd i)D R d R + h.c, (5) 



TABLE II: Particle contents for the model F 



Fields 


SU(3) C 


SU(2) L 


U(1) Y 


U(1)g 


Ql 
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+1/6 


+1 


Ur 


3 
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+2/3 


+1 


Dr 


3 


1 


-1/3 


+1 


il> 


1 


1 
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where H = ia 2 H*, Q L = (U L D L ) T , q L = (ulcI l ) t is the SM doublet quark of the 1st 
generation, and Au,d are parameters with mass-dimension one. 

After the electroweak symmetry breaking, the mass eigenstates of U are obtained as 

U L \ = I cos9 u sin9 u \ fuA 
U R J y- sin 9 U cos 9 U J \U 2 J 

with a mixing angle 9 U . Similarly, D\ and D 2 are obtained with a angle 9 d . With the mass 
eigenstates, the Yukawa interactions between the dark matter fermion and the SM quarks 
in Eq. (jSJ) are rewritten as 

£ Y = -Tp(f L cos 9 U P L - f Ru sin B u Pr)U\u - 4>{f L sin 9 U P L + f Ru cos 9 u P R )U\u 

- $(f L cos6 d P L - f Rd sm9 d P R )D\d-$(f L sm9 d P L + f Rd cos9 d P R )tyd + h.c. (7) 



III. DARK MATTER ELASTIC SCATTERING WITH NUCLEI 



The dark matter scattering cross section with nucleus (N) made of Z protons (p) and 
A — Z neutrons (n) is given by 

1 / rriN 



a. 



(Zf p +(A-Z)f n ) 



TV \m N + m^ 
for a scalar dark matter, while for a Dirac fermion dark matter 

2 



a. 



1 / m N m^ 



(Zf p + (A-Z)f n ) 



(8) 



(9) 



The effective coupling with a proton f p and a neutron f n is expressed, by use of the hadronic 
matrix element, as 



m, ^ JT9 m 27 iTG ^m 

q=u,d,s H c,b,t H 



(10) 



where a q is an effective coupling of the DM particle with a g-flavor quark defined in the 

operators 

J a q qq\(j)\ 2 for 

^-int — < , (.J- J- J 

I a q qqipip f° r ^ 
with its mass m g , /^V and /^ with i — p,n are constants. In our analysis, we use the 
following values: /® = 0.0269, /« = 0.0390, /g = 0, and /$, = 1 - E«=„,d, 8 /$■ Those 
/y„ and /^ are quoted from Ref. [34| . while we set fj^ s = because recent studies of the 



lattice simulation 35] as well as chiral perturbation theory 



imply negligible strange 



quark content. It has been pointed out |Tj| that the results of XENON100, CoGeNT and 
CRESST can be compatible, if the following relations are satisfied: 

j- ~ -0.7, ag,-2x 10~ 2 pb. (12) 

Jp 

Note that /„ ^ f p , and therefore the dark matter particle has isospin violating interactions 
with quarks. 

1. Model S 

For the model S, there are two contributions to the effective coupling a q . One is from 

the exchange of the scalars h\ and hi, for which we find 

A sin 2a / 1 1 \ A sin 2a 

a q = -m q — r ~ -m q — — — , (13) 

where we have assumed raL <C m \ 2 , with mh 2 being the SM(-like) Higgs boson mass. Note 

that a q /m q is independent of q, so that this contribution conserves the isospin. The other 

contribution is from the exchange of the Dirac fermions, U and D: 

fj m fl m^ ■ ^ fhrrn fhrn^d 

q 2 M 2 -m\ q + 2 M 2 D -m\ q ~ 2 M 2 q + 2 M 2 D q ' [ } 

where we have assumed m\ <C My D . Clearly this contribution violates the isospin symmetry. 
For simplicity, let us assume fp -C fu, and the total contribution is given by 

««- rn q 2m ^ + 2 ^d,. (15) 

Note that the existence of the two terms (with A sin 2a > 0) is crucial to realize the op- 
posite signs between f p and f n , because the heavy quark U (and also D) always positively 
contributes to a„. 
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FIG. 1: The contours of scattering cross section with a proton for various values, o"gj = 0.02 pb 
(solid), 0.1 pb (dashed) and 0.5 pb (dotted), together with the (red) straight line along which the 
condition f n /f p ~ —0.7 is satisfied. Here we have fixed the dark matter mass as m^ = 8 GeV. 



Fig. Q] shows the contours for various values of Cgj, along with the (red) straight line 
corresponding to the condition f n /f p = —0.7. The two conditions in Eq. ( 1T2|) are satisfied 
for 

\/A sin a cos a „ „ A „^_ a ^ Tr _, |/; 



3.14 x 10" 2 GeV" 1 , 



u\ 



rn hl M v 

Here we have fixed the dark matter mass as irtj, = 8 GeV. 



3.03 x 10" d GeV 



r-l 



(16) 



2. Model F 



For the model F, the effective coupling a q is given by 



a„ 



sm26 u f L f Ru 



M 2 



M 2 f "'' 
u 2 



6^ + sm29 d f L f Rd 



M k 



M k 



sm26 u f L f R sin26 d f L f R d 

Ml q+ M 2 q 



(h 



Di 



(17) 



Here, for simplicity, we have taken a limit, M 2 ~ -C M 2 ~ and M 2 ~ <C M 2 ~ . This effective 

1 J1 Ui U2 Di D2 

coupling violates the isospin symmetry and f n / f p < can be realized when the relative signs 
between sin 29 u f Ru and sin 29 d f Rd are opposite. We further simplify the system by setting 



f L cos 6 U = f Ru sin 6 U = ffj > 0, 
f L cos6 d = -f Rd sin6 d = f D > 0, 



(18) 



so that 



a q ~ 



fu \ 5 u + (jj^} dt 



1V1 Ui/ \ IV1 D 1 



(19) 



Fig. |2] shows the contours for various values of <7gj, along with the (red) straight line 
corresponding to f n /f p = —0.7. The two conditions in Eq. ( !T2|) are satisfied for 



/, 



u 



M, 



1.58 x 10" d GeV 



/. 



D 



U x 



M 



1.82 x 10 _d GeV 



(20) 
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IV. DARK MATTER ANNIHILATION CROSS SECTION 



In this section, we estimate the annihilation cross section of the scalar /fermion dark 
matter particles for the parameters identified in the previous section to satisfy the conditions 
for the IVDM. We will see that the s-wave annihilation cross section of the scalar dark matter 
is too large to reproduce the observed relic abundance while that of the fermion dark matter 
appears to be of the right order of the magnitude. However, both of those models cannot be 
viable, because the s-wave annihilation cross section exceeds the upper bound obtained by 
the Fermi-LAT observations [16| . In the last part of this section, we will propose a simple 
scenario to realize the IVDM being consistent with the Fermi-LAT observations. 
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FIG. 2: The same as Fig. Q] but for the Dirac fermion dark matter. Here Xjj 
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A. Model S 



The dominant dark matter annihilation process is found to be <f>'<f> — > bb mediated by the 
scalars, hi and h 2 , in the s-channel. Assuming vr? hx < m 2 , <C m2 h 2 , the s-wave annihilation 
cross section is evaluated as 

. . 3 / A sin 2a m b \ 

( "> = w, { i^^i) • (21) 

where mj = 4.2 GeV is the bottom quark mass. Using the values in Eq. f fT6|) . we find, 
for example, (av) ~ 6.5 pb for m^ = 8 GeV and m^ = 3.9 GeV. This cross section 
is roughly one order of magnitude larger than the typical dark matter annihilation cross 
section (av) ~ 1 pb to achieve the observed relic density. Thus, in this case, the resultant 



dark matter abundance becomes too small. In order to realize the observed relic density, 
we may assume a nonthermal production of dark matter particles in the early Universe. 
However, this cannot be a phenomenologically viable scenario, because the dark matter 
annihilation cross section to the bottom quarks is constrained by the Fermi-LAT data as 
(av) < 0.5 pb [la]. Practically, this constraint can be satisfied if we take m^ < 2.1 GeV. In 
this case, the dark matter is overabundant and the relic density should be diluted by some 
mechanism in the history of the Universe. Since such a scenario is quite ambiguous, we do 
not consider it in this paper. 

B. Model F 

The s-wave annihilation modes are given by t-channel U/D exchange with uu/dd final 
states. In a limit m 2 <C M% - , the cross sections is found to be 

V Ul,2,i->1,2 



(av) ~ a u + a d , 



(22) 



with 
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An 
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H TT75 I + Jr. I ^TTT5 r 
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(23) 



(24) 



where we have used Eq. (fTHD and the limit M% < M% and M? <C M? . Using the values 

y ' Ui U2 i/1 i<2 

in Eqs. ( 120]) . we find the annihilation cross section as 



(av) ~ 0.61 pb. 



(25) 



Similarly to the Model S, this annihilation cross section is excluded by the Fermi-LAT 
observations 16(. The upper bound on the cosmic anti-proton flux obtained by the Fermi- 
LAT observations is interpreted to a cross section upper bound of DM annihilations to up 



10 



and down quarks as [26 1 



{av) < 0.2 pb. (26) 

C. Solution to too large annihilation cross section 

As we have seen, for a given parameter set to realize a large enough isospin violating 
scattering cross section with nuclei, the resultant annihilation cross section is too large 
to satisfy cosmological and astrophysical constraints. For relic density, one may assume 
a nonthermal dark matter production. However, as we have seen, such an idea cannot 
work because of the severe upper bound on the dark matter annihilation cross section from 
the Fermi-LAT observations. In order to avoid the Fermi-LAT constraints, we propose an 



extension of our model to the so-called "asymmetric dark matter" scenario 37H43|. This 
scenario is suitable to our model, because the global U(l)c symmetry introduced in our 
model leads to the conservation of the dark matter number. Once a suitable DM-anti- 
DM asymmetry is created in the early Universe, the too large annihilation cross section 
nicely works to leave only the dark matter in the present Universe with the observed relic 
abundance. Since the relic abundance of anti-dark matter particles in the present Universe 
is much smaller than the dark matter one, a cosmic ray flux produced by DM and anti- 
DM annihilations becomes much smaller and hence the constraint from the Fermi-LAT 
observations can be avoided. 

A relic density of the dark matter particles in the presence of dark matter asymmetry 



(chemical potential) has been analyzed in detail by solving the Boltzmann equations 45[. For 
example, with a suitable initial dark matter asymmetry, the observed relic abundance of the 
dark matter particle can be obtained by the s-wave annihilation cross section (av) = 0(1) 
pb, while the relic abundance of anti-dark matter particle is found to be two orders of 
magnitude smaller than the dark matter one. As annihilation cross sections become larger, 
the relic abundance of anti-DM particle becomes exponentially smaller. This result is almost 
independent of WIMP dark matter mass. We apply the result to our scenario, so that the 
cosmic ray flux from DM-anti-DM pair annihilations is significantly suppressed and the 
constraint from the Fermi-LAT observations is avoided. 
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V. CONSTRAINTS FROM COLLIDER EXPERIMENTS 

A. Constraints on the mediator (s)quarks from LHC 

Our model includes heavy (s)quarks, which can be produced at the Large Hadron Collider 
(LHC) mainly through the gluon fusion process. The heavy (s)quarks, once produced, decay 
to the SM quarks and the dark matter particles, and this process is observed as the hadronic 
final states with transverse missing energy. Searches for such events have been performed 
at the LHC experiments, and the null result, so far, sets the lower bound on heavy (s) quark 
masses as > 800 GeV |44j. This bound is obtained for the so-called simplified MSSM, 
where scalar quarks of the first two generations are produced at the LHC and decay to 
quarks and the lightest superpartner neutralino. Since we only introduced one generation 
of heavy (s) quarks, the mass bound on the mediator (s) quarks should be a little milder, 
but let us apply the bound for conservative discussion. From Eqs. ( TT6|) and (|20|) . we can 
see that this mass bound is satisfied with the couplings being in a perturbative regime, 
^/(4vr),/?y(4vr)«l. 

B. Constraint from Higgs boson invisible decay 

In Model S, the scalar mass eigenstate h 2 is approximately identified as the SM Higgs 
boson. Through the mass mixing with the singlet scalar s, the SM Higgs boson decays to a 
pair of the dark matter particles 1 . This decay width is given by 



X 2 v 2 sin 2 a / 4m 



The current ATLAS 1471 and CMS 



^"^-nsssrV 1 <' (27) 



481 ] data for the Higgs boson production and its various 



decay modes are mostly consistent with the SM expectations, and the branching ratio of 



invisibly decaying Higgs boson is constrained (at 3a) as 49] 



BR(/i 2 ->■ invisible) = l~* — ' < 0.35, (28) 

Tsm + r(/i 2 — > 



where r SM = 4.07 MeV [50j is the SM prediction of the total decay width of a Higgs boson 
with a 125 GeV mass. 



1 This structure is the same as in the so-called Higgs portal dark matter scenario. See, for example, Ref. [4; 
for detailed analysis and references therein. 
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Using the result in Eq. ( 1T6|) . we can give the annihilation cross section of Eq. (12 ip and 
the Higgs invisible decay rate of Eq. (128]) as a function of only m^ , with a fixed dark matter 
mass m<£ = 8 GeV. The correlation between these two quantities is shown in Fig. [3] by 
varying rrih 1 in the range of 1.4 GeV < m^ < 7.0 GeV. Here the vertical line denotes the 
upper bound, BR(/i 2 — > invisible) = 0.35 at 3cr [49j while the horizontal line corresponds 
to a typical value ((crv) = 1 pb) of the WIMP dark matter annihilation cross section for 
reproducing the observed relic abundance. The upper bound BR(/i 2 — > invisible) = 0.35 is 
obtained by m^ ~ 3.9 GeV, for which we find the annihilation cross section (erf) ~ 6.5 
pb. Note that the asymmetric IVDM scenario we have proposed in the previous section 
can be consistent with the constraint on the Higgs invisible decay rate. In order for the 
asymmetric dark matter to be consistent with the observed relic abundance, we have a 
lower bound on the annihilation cross section as (crv) > 1 pb |45J. Applying this bound, we 
read BR(/i 2 — > invisible) > 8% from Fig. [31 Precision measurements of Higgs decay width 
at future collider experiments such as the international linear collider, photon collider and 
muon collider can reveal the existence of the dark matter. 

VI. CONCLUSIONS 

The IVDM scenario with destructive interference between the dark matter scatterings 
with proton and neutron offers an interesting possibility to reconcile conflicting results among 
direct dark matter search experiments for a light WIMP with mass around 10 GeV. In this 
paper, we have considered two simple IVDM models and investigate various phenomeno- 
logical aspects of the models, such as realization of the IVDM scenario, the constraints on 
dark matter annihilation cross sections from the dark matter relic abundance as well as 
indirect search for dark matter, and collider experimental constraints on the extra particles 
introduced in our models. 

One model introduces a complex scalar as dark matter particle along with heavy extra 
quarks and a SM singlet real scalar, through which the dark matter particle couples with the 
SM up and down quarks. Isospin violating effective couplings are realized by the interference 
between processes mediated by the heavy quarks and the scalar. In the other model, we 
have introduced a Dirac fermion as dark matter particle along with heavy colored scalars 
analogous to squarks in the MSSM, through which the dark matter particle couples with 
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FIG. 3: The correlation between the Higgs invisible decay rate and the dark matter annihilation 
cross section through 77i/ ll . Here we have varied m/ ll in the range of 1.4 < m/ ll (GeV) < 7.0. 

the SM quarks. The interference between two processes mediated by up-type squarks and 
down-type squarks realizes the isospin violating effective couplings. For both models, we 
have identified a parameter region suitable for the IVDM scenario. With the parameter 
regions, we have also calculated the dark matter annihilation cross section to examine the 
thermal relic density and consistency with indirect dark matter searches. 

We have noticed that for both models, the calculated dark matter annihilation cross 
sections exceed the upper bound obtained by the Fermi-LAT observations, and therefore 
the parameter regions for realizing the IVDM scenario are excluded. We have proposed a 
simple scenario to reconcile the IVDM scenario with the Fermi-LAT observations, namely, 
the IVDM being asymmetric. In our models, a global U(1)q symmetry has been introduced 
whose conservation ensures the stability of dark matter particle. At the same time, this 
global symmetry leads to the conservation of the dark matter number and this structure 
is suitable for the asymmetric dark matter scenario. As discussed above, we have found 
that the dark matter annihilation cross section is too large to satisfy cosmological and 
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astrophysical constraints simultaneously. In fact, when a suitable asymmetry between dark 
matter-anti-dark matter is generated in the early Universe, the large cross section nicely 
works to leave only the dark matter in the present Universe. Thus, the relic abundance of 
anti-dark matter particle is much less than the dark matter relic abundance, as a result, 
the flux of cosmic rays created by annihilations of the dark matter and anti-dark matter 
particles is dramatically suppressed and the constraint by the Fermi-LAT observations is 
avoided. Since a variety of models to account for generating the dark matter asymmetry 



has been proposed (for an incomplete list, see e.g., [5lJ-[60J]), we do not propose a specific 
mechanism for it in this paper. 

We have also considered collider experimental constraints on our model. Colored fermions 
and scalars introduced in our models can be produced at the LHC and their decays to 
the SM quarks and dark matter particles yield the signal events with jets and missing 
transverse energy. We have confirmed that our IVDM scenario is realized consistently with 
the current LHC bound on the mass of the colored particles. In the Model S, the SM 
Higgs boson invisibly decays to a pair of dark matter particles and the upper bound on the 
invisible decay rate is given by the LHC data. We have found a parameter region in which 
the IVDM scenario is consistent with the LHC bound on the Higgs boson invisible decay 
rate. Interestingly, our successful asymmetric IVDM scenario leads to a lower bound on the 
invisible decay rate about 8%, so that precision measurements of the Higgs decay width at 
future collider experiments can test our scenario. 

Observable effects of the asymmetric dark matter scenario in neutron stars have been 



investigated [61 



621 ] . Since the dark matter particles do not self-annihilate, once captured 
in neutron stars, dark matter particles are continuously accumulating and neutron stars 
eventually collapse into black holes. Observations of old neutron stars provide constraints 
on parameters of the asymmetric dark matter scenario. In particular, such constraints are 
more severe for the case with a scalar dark matter because of the absence of Fermi degeneracy 
pressure. However, since the resultant constraints highly depend on the strength of dark 
matter self-interactions [62], we do not consider the constraints from the black formation in 
our scenario. 
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